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Abstract 

We study the spectral statistics in the center of the lowest Landau band of a 2D disordered system with smooth 
potential and strong transverse magnetic field. Due to the finite size of the system, the energy range in which there 
are extended states is finite as well. The behavior in this range can be viewed as the analogue of the Anderson metal- 
insulator transition for the case of the Hall system. Accordingly, we verify recent predictions regarding the exponent 
of the asymptotic power law of S 2 (7V), 7, and that of the stretched exponential dominating the large s behavior of 
the spacings distribution, a. Both the relations, a = 1 — 7, and j = 1 — where v is the critical exponent of the 
localization length and d is the dimension, are found to hold within the accuracy of our computations. However, we 
find that none of several possible models of the entire spacings distribution correctly describes our situation. Finally, 
for very large N, N > 60, we find a new regime in which S 2 (7V) behaves as a power law with an unexpectedly large 
power, 71 = 1.38 ±0.02. 
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1. Introduction 

The microscopic description of disordered systems often 
relies on various types of Random Matrices. These are 
in general different from the traditional Gaussian Ran- 
dom Matrix Ensembles, GOE, GUE and GSE, which 
were introduced in the context of nuclear physics in the 
early fifties. The Gaussian Ensembles are mainly used 
to describe the statistical properties of the spectra of 
strongly interacting, complex quantum systems such as 
nuclei for example. Consequently, such ensembles are as- 
sumed to be invariant under similarity transformations 
and therefore, their members are basically structureless 
matrices with elements, hij , that are independent, Gaus- 
sian distributed random variables. On the other hand, 
in disordered systems, the central phenomenon is the lo- 
calization of eigenvectors in configuration space. While 
such phenomenon is in general absent in the Gaussian 
Ensembles, it is characteristic of random matrices with 
off-diagonal elements that decay as a function of \i — j\. 
Such matrices have a preferred representation and are 
often referred to as banded. 



In the definition of the Gaussian Ensembles all the 
structure of the physical system is ignored, except for its 
symmetry with respect to time reversal. Accordingly, 
their prediction power is based on the assumption of 
universality, although the extent of the corresponding 
universality class is only qualitatively determined. For 
example, disordered systems in the metallic regime have 
been shown to share most of the spectral properties of 
the Gaussian Ensembles for small enough energy inter- 
vals. In particular, the spacings distribution, P(s), is 
of Wigner type and the number variance, £ 2 (7V), is log- 
arithmic. Disordered systems in the insulating regime 
however, belong to a different universality class. Namely, 
their P(s) is Poisson and Y, 2 (N) = N for N < N T . In 
order to understand the spectral characteristics of the 
insulating regime, one can think of the sample as being 
composed of subsystems of the size of the localization 
length, £. While the spectrum of each subsystem is of 
metallic type, the full spectrum consists of a random 
superposition of many such spectra and therefore is en- 
tirely uncorrelated, or in other words, Poissonian. The 
universal behavior in both the metallic and the insulat- 
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ing regimes is restricted to a certain, system dependent, 
energy range known as the Thouless energy, Et = Ti/t, 
where r is the time it takes an electron to diffuse through 
the sample. 

It was recently suggested that a third type of spec- 
tral behavior which is neither Wigner nor Poisson should 
be expected in the neighborhood of a metal-insulator 
(MI) transition. 1,2 Using a perturbative approach, it was 
shown that for a rf-dimensional system, inside an energy 
interval, 6E, which is centered on the transition energy, 
E c , and in which £ > L, where L is the system size, 

P(s) oc e^i-AdPs 2 " 1 ) for s -> oo , (1) 

and 

£ 2 (#) oc for N oo . (2) 

P 

Iff is the critical exponent of £, £ cx (E—E c )~ v , then 7 = 
1 — Finally, j3 = 1,2,4 according to the symmetry 
with respect to time reversal. 

Several numerical studies of the spectral statistics in 
the neighborhood of the MI transition in the 3D Ander- 
son model were also performed. In particular, in Ref. 3 
the family of P(s) functions that interpolate between the 
Wigner and the Poisson form is obtained for a sample of 
finite size. All the curves in this family cross in one point 
at s ~ 2 and their maxima appear to be equally high. 4 
Moreover, P(s) = c\s for s<l and for the critical distri- 
bution (W = 16), c\ ~ 2.04. Scaling theory implies that 
the family of P(s) is parametrized by a single variable, 
P(s, L/Zoo), where ^ = lim i _ 00 £ i . In the thermody- 
namic limit, L — ► 00, only the Wigner, the Poisson and 
the critical P(s) can occur. As the size of the system 
grows, the other curves gradually migrate towards the 
Wigner distribution for W < W c and toward the Pois- 
son P(s) when W > W c . In another study, Evangelou 5 
has shown that 

P v (,s) = Cl /exp(- C2S 1+a ) , (3) 

with (3=1 gives the best fit to the numerically obtained 
critical P(s) for c\ = 2.65, c 2 = 1.47 and a = 0.31. 
Moreover, the fitted curve appears to be in good agree- 
ment with the numerical one. For the value of a the 
error was estimated around ±0.06. Consequently, it was 
suggested that for v = 1 this confirms the prediction of 
Eq. (1), a = 0.33. Although the value oiv is not yet fully 
agreed upon, most studies indicate that v ~ 1.35. This, 
in turn, leads to a ~ 0.25 which still agrees with the nu- 
merical value. When instead a power law was fitted to 
In P(s) for 2 < s < 4, the resulting a was quite low, that 
is, a ~ 0.12. The behavior of the critical £ 2 (7V) is also 
shown in Ref. 5. It indeed displays power law behavior 



with 7 ~ 0.88. It therefore appears that there is better 
support for the relation 

7 = 1 - a , (4) 

where the large s value of a is considered than for the 
relation between 7 and v. As was pointed out in Ref. 6, 
Eq. (4) holds quite generally. In a large enough energy 
interval, s, the number of levels, N , can be assumed to 
be Gaussian distributed. If the corresponding average is 
N and the variance is S 2 then the probability of having 
no levels in s, P(s), is 

P(s) oc exp(-# 2 /S 2 ) = exp(-cs 2 - 7 ) , (5) 

where Eq. (2) was used for S 2 and s is in units of the 
mean spacing. 7 Therefore, Eq. (4) also holds for separate 
ranges of N and the corresponding ranges of s. Another 
feature of Y, 2 (N), shown in Fig. 2(b) of Ref. 5 but oth- 
erwise ignored, is that a power law behavior is observed 
also for very small N although with a slightly smaller 
exponent, 70 ~ 0.71. The change in the exponent occurs 
at N = No ~ 2.5. Moreover, it is interesting to notice 
that 70 and the value of a obtained from fitting Eq. (3) 
to the data satisfy Eq. (4) well within the correspond- 
ing errors. As we shall see in the next Section, such fit 
is dominated by the s < 2 range and accordingly, the 
resulting a is determined by the same range of N as 70. 
On the other hand, it is surprising that the Gaussian ap- 
proximation should hold here. Another fit of the critical 
P(s) with Eq. (3) was done in Ref. 8. Since here the 
analytical forms of c\ and c 2 obtained from the normal- 
ization and s = 1 requirements were used, only a one 
parameter fit was performed. The resulting exponent 
was a = 0.20 ± 0.03 which corresponds to a slightly too 
large value of v , v = 1.7 ± 0.2. It is not clear whether 
the difference between the value of a obtained in Ref. 5 
and that of Ref. 9 is due to the quality of the fit or the 
different parameters of the model. 

In this paper we study the level statistics in the neigh- 
borhood of another MI transition, namely, that occur- 
ring in two dimensional disordered systems in the pres- 
ence of a strong transverse magnetic field. Such systems 
have been found to display quantized Hall conductance 
which is a consequence of the divergence of £ at the cen- 
ter of the disorder broadened Landau level. There is 
considerable evidence that the corresponding critical ex- 
ponent is v = I and that this value is more reliable than 
the one for the Anderson transition. Therefore, the Hall 
system is preferable for verifying the predictions of Eqs. 
(1 - 2). On the other hand, since here extended states 
can only be found precisely at the center of the band, it 
is not clear whether these predictions apply. However, 
in a system of finite size the critical region is broadened 
into a range of energy, AE, where £ > L. Inside AE 
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the Hall system is critical in a way similar to that of the 
Anderson type system inside 6E. 

Numerical studies of the corresponding spectral statis- 
tics of Hall systems 9 were done prior to the work of 
Refs. 1 and 2. It was found that the spacings distri- 
bution varies from an almost GUE P(s) in the center 
of the Landau band to a Poisson one at its edges. The 
corresponding transition was also observed in the behav- 
ior of the spectral rigidity, Aa(N). On the quantitative 
side, it was suggested that the ensemble of Pandey and 
Mehta 12,13 that interpolates between GOE and GUE can 
reproduce the P(s) obtained close to the center of the 
band. The members of this ensemble, Hpm, are of the 
form 



Hpm(oipm) — Hs + iapAiHA, 



(6) 



where Hs belongs to the GOE matrix, Ha is a real anti- 
symmetric random matrix and < «pm < 1 such that, 
Hpm{X) = GUE. Presently, only the P(s) for the case 
of 2 x 2 matrices is known 



PM 




(7) 



where v is determined by the requirement that s = 1, 




apM 



: arctan 
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l PM 



apM 
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In particular, for the P(s) calculated from the first 20 
levels around the center of the band (out of a total of 200 
levels in each sample and with 700 samples being used) 
good agreement was obtained with Ppm{®-pm = 0.515). 
Notice however that for large s the error function in Ppu 
approaches unity and thus j = 0, in contradiction with 
the prediction of Refs. 1-2. Qualitative studies of the 
P(s) for Hall systems were also reported in Refs. 10 and 
11. 



2. Results 

In order to obtain the spectrum of an electron in a 
2D random potential and a strong magnetic field, B = 
Bz, we numerically diagonalize the Hamiltonian matrix, 
Hkk 1 , in a lowest Landau level basis, <S>j.(x,y). In the 
Landau gauge, 



(9) 



where lp = ^hc/eB is the magnetic length and 
v ( x > y) = ^2 v nm exp ^— — X "* ^J; V — 



(10) 



The v nm coefficients are uncorrelated random numbers 
chosen from a uniform distribution with zero mean and 
variance u> 2 /12 and (x n , y m ) = (— a/2 + (n — I) a, —a/2 + 
(m — l)a) where a is the lattice constant. Moreover, 
the system is restricted to a square of size L such that 
Ni = L/a is an integer and accordingly, 1 < n,m < 
The corresponding Lowest Landau Level (LLL) states 
are A'-fold degenerate in the absence of disorder, K = 
L 2 /(2iTlg), and are chosen to satisfy periodic boundary 
conditions 14 



<t> k (x,y) = (Li B r 1/2 *- 1/4 £ 



k/L+jL/ll)x] 



x e 



-[(y-jL)/l B -2Trkl B /L] 2 /2 



(11) 



where k is integer and —K/2 < k < K/2. Finally, 
shifting the origin of the energy such that the LLL is 
at E = 0, one obtains 



H. 



1 



kk' 



(12) 



where 



Sw = -^B + 2<7 2 )fe - q' k ) 2 + i(q k - q'k)x r , 



ik + <i' k ,2 

7. ' R 



(13) 



and qk = ^f-k. Since the terms in the second sum of 
Eq. (12) decay rapidly with j and j' , it is sufficient to 
truncate the sum such that —1 < j,j' < 1. The parame- 
ters determining the strength and shape of the disorder 
are w, a and a. We first require that the overlap be- 
tween two adjacent impurities be exp(— 77). Specifically, 
consider the overlap integral 



R[a,a)= j dxdyexp 



x 2 + y 2 



exp 



2a 2 
(x - a) 2 + y 2 
2a 2 



(14) 



Then 



R(a, a)/R(0, a) = exp ( = exp (-^ 



(15) 



3 



which leads to 



(16) 



Moreover, we assume a = lp such that there is one impu- 
rity within an area of l B . Let us now fix the strength of 
the disorder, w. Notice that after projecting the Hamil- 
tonian on the lowest Landau level, the value of w only 
determines the energy scale of the problem. For conve- 
nience we let the width of the density of states, ^2, be 
unity 



/ u 2 (a 2 = 2a) 



24 7 ra 2 (/ 2 f + a 2 ) 



1 



(17) 



Finally, the value of B in our computations corresponds 
to K = 1019 and we have used 160 different samples. In 
order to avoid the tails of the Landau band, we restrict 
ourselves to the energy range |-E"| < E m =2.5, which in 
turn we divide into Np = 5 equal intervals. In the third 
interval, \E\ < 0.5, both the density of states and the 
localization length are almost constant. Therefore, in 
what follows, we shall regard the properties of the spec- 
trum in this interval as critical. In order to compute the 
spacings distribution, we have first unfolded the spec- 
trum normalizing each individual spacing to the average 
of the Np = 5 adjacent ones. In Fig. 1 we show the 
resulting P(s) for the case where the Np = 119 spacings 
closest to E = were used. Three different theoreti- 
cal curves are fitted to the numerical experiment: (a) 
the distribution of Eq. 3 with [3 = 2, Py(s), (b) the 
Pandey-Mehta distribution of Eqs. (7 - 8), Ppm(s), and 
(c) the Robnik distribution 15 , Pr(s). 

The Robnik distribution interpolates between the 
Poisson P(s) and the GUE one using a 2 x 2 random 
matrix, Hr(c(r). For cur = 0, the matrix is diagonal and 
P(s) = e~ s , while for cur — ► 00 it approaches the GUE. 
In practice, when cur ps 1 the corresponding spacings 
distribution is very close to that of GUE. Specifically, 



P R (s) = P+(s) + P-(s), 
where 

8 e X2 r S ' a P T((S/a) 2 -^ 2 ) 1/2 

P±(-s) = -,- T 



(18) 



-j; 2 /4A 2 



4a 2 A 2 J ((S/a) 2 - x 2 ) 1 ! 2 
x [1 ± erf (^((S/a) 2 - £ 2 ) 1/2 T \)]dx, (19) 

A = urI a and a is determined by the s = 1 constraint. 

Fitting Pv(s) to the numerically obtained P(s) with 
the assumption that the number of spacings in each of 
the 38 bins is a Gaussian distributed random variable 
leads to a = 0.26 ± 0.02 and X 2 = 85.0 (see Fig. 1(a)). 



For the case of a least squares fit with 37 degrees of free- 
dom this corresponds to a negligible confidence limit, 
CL = 0.7 • 10 _3 %. One might wonder whether the large 
value of x 2 obtained is a consequence of ignoring addi- 
tional sources of error besides those of statistical origin. 
The most likely such source is the unfolding procedure. 
In order to verify this possibility we modify the extent 
of the smoothing interval, Np, to 7. Comparing the re- 
sulting P(s), fV(s), against the original one, -Ps(s), and 
using twice the statistical variance we obtain \ 2 = 29.1 
for which CL ~ 62%. Therefore, despite significantly 
larger statistics than usually encountered in Quantum 
Chaos where in general only a single sample is studied, 
here as well, the accuracy of the unfolding procedure ap- 
pears to be practically irrelevant. One concludes that 
Pv(s) is inappropriate for the description of our numer- 
ical results. Nevertheless, it is interesting to notice that 
most of the contribution to the \ 2 comes from the small 
s range. In particular, we find that the partial \ 2 from 
the first 10 bins, s < 1, is 48.2 while that from the last 
18 bins, s > 2, is only 18.0. This fact suggests that the 
large s tail of Pv(s) could be in agreement with P(s) 
despite the fact that the two distributions are different 
from each other. We therefore separately fit the s > 2 
tail of P v (s) and this leads to a = 0.30 ±0.02, X 2 = 13.9 
and CL ~ 74%. This result is statistically reliable and, 
at the same time, is in contradiction with the prediction 
of Ref. 2 for v = 7/3 which is a = 0.214. Alternatively, 
using this result to determine the value of v one obtains, 
v = 1.7 ± 0.1. 

While Pv(s) does incorporate the expected behavior 
of the true P(s) at both small and large s, it is largely 
an interpolation formula for s = 0(1). In other words, 
there is no physical argument to support such behav- 
ior as that of Pv(s) for intermediate values of s. On 
the other hand, the Pandey-Mehta distribution interpo- 
lates between the GOE P(s) and the GUE one along the 
most natural path in the space of 2 x 2 random matrices. 
At first, it would appear that there should be no GOE 
component in the description of a Hall system. How- 
ever, such component is needed in order to incorporate 
those states that are extended, £ > L, but are relatively 
insensitive to changes of the flux and accordingly, be- 
have almost as if time-reversal invariant. We therefore 
fit Ppm(s) to the numerical P(s) (see Fig. 1(b)) and ob- 
tain apra = 0.18 ± 0.01 and \ 2 = 94.4. In other words, 
the Pandey-Mehta distribution is even less appropriate 
that Pv(s) to describe the level spacings at the center 
of the Landau band. Moreover, we find that the contri- 
bution to the large value of \ 2 comes equally from the 
entire range of s. Specifically, the partial \ 2 f° r small s, 
s < 1 (the first 10 bins), is 23.3, while that from large s, 
s > 2 (last 18 bins) is 54.0. Accordingly, when the fit is 
restricted to the large s tail, s > 2, we obtain otpM = 
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and x 2 = 37.2. Not only does this correspond to a con- 
fidence limit of only 0.32%, but rather than a true min- 
imum of x 2 it is just the edge of the allowed interval of 
the fitting parameter, «pm- One is lead to conclude that 
the tail of the distribution is not of Gaussian type, that 
is, a ^ 1. 

If instead of computing the spacing distribution in the 
central (third) energy interval, we use the fifth, edge in- 
terval, a P(s) that is quite close to a Poisson distribu- 
tion is obtained (see Fig. 2). This is to be expected 
considering that the corresponding eingenstates are lo- 
calized, ^ <C L. On the other hand, from the knowledge 
of the finite size localization length, ^r, alone, there is 
no way one can determine the boundaries of the criti- 
cal energy range. It is therefore conceivable that even 
in the central interval some fraction of the states are 
localized. Since the Pandey-Mehta distribution ignores 
the possibility of a Poisson component due to localized 
states, it is worthwhile to try to identify the presence of 
such component by fitting a Poisson - GUE interpolat- 
ing distribution, e.g. Pr(s) (see Fig. 1(c)). We obtain 
that cxr = 0.66 ±0.02 and x 2 = 604.4 which corresponds 
to a negligible confidence limit. Although here, like for 
Py(s), most of the contribution to x 2 comes from the 
small s range (x 2 (s < 1) = 365.2, X 2 (s > 2) = 106.8), 
the fit of the s > 2 tail alone leads to cur = 0.75 ± 0.02, 
X 2 = 92.8 and CL < 10" 3 %. This clearly excludes the 
Robnik distribution as a candidate for the description of 
our results. Moreover, it indicates that in order to derive 
an appropriate distribution one has to use the fact that 
some of the eingenvectors are localized. It is likely that 
such distribution is similar to the one encountered in a 
ID disordered system when ^ L (see Refs. 16 - 20). 
The failure of Pr(s) to incorporate localization is further 
stressed by fitting it to the edge P(s) of Fig. 2. The best 
fit is obtained for cxr = 0.75 • 10 -3 with x 2 = 85.9 and 
CL < 10" 3 %. 

As already mentioned, we do not have a clear estimate 
of the critical energy range. It is therefore instructive 
to verify the behavior of the tail of P(s) as we shrink 
the energy interval around the center of the band. We 
therefore progressively reduce the number of energy lev- 
els around E = 0, Nr, and fit Py(s) for s > 2. The 
resulting values of a are given in Tbl. 1. One would 
expect that as we lower Nr and correspondingly, the 
statistics becomes worse, the value of x 2 would go down 
and the error would increase. We see however that while 
the error does indeed increase, the x 2 goes first way up 
before returning to about the same value at Nr = 9 
as at Nr = 119. Therefore, although in Tbl. 1 the 
value of a is varying beyond what our error estimates 
would allow, it is reasonable to expect that the values 
at Nr = 9 and Nr = 119 are the most reliable. Using 
the former together with the prediction of Ref. 2 to ob- 



tain the critical exponent, leads to v = 1.7 ± 0.4, which 
still disagrees with the corresponding theoretical value, 
v = 7/3. However, we incline to view this disagreement 
as being a consequence of our inaccurate numerical pro- 
cedure rather than indicating the failure of the theory of 
Ref. 2. 

We now turn to the study of £ 2 (7V). In Fig. 3 we 
show the result obtained for the third, central interval. 
Ignoring the range, 0.8 < N < 2.3 in which some oscil- 
lations occur, one can distinguish three different power 
law ranges: I. 0.07 < N < 0.8, II. 2.3 <N < 60, and III. 
60 < N < 350. The best fitting power law curves give 
70 = 0.74 ± 0.01, 7 = 0.77 ± 0.02 and 7l = 1.38 ± 0.02, 
respectively. Moreover, the corresponding natural loga- 
rithms of the coefficients are Co = —0.69 ± 0.02, C = 
-1.05 ± 0.04, and C\ = -3.5 ± 0.1, respectively. As in 
Ref. 5, 70 and the value of a obtained from the fit of 
the entire Py satisfy Eq. (4). This relation also holds 
for 7 and the a of the s > 2 tail of Py but only if the 
large error bar of the Nr = 9 case is used. On the other 
hand, the relation between 7 and v is in good agreement 
with the prediction of Ref. 1. The value of 71 however, 
bears practically no relation to any previous result or 
theoretical prediction. While presently the origin of this 
regime is not clear, we have observed similar behavior in 
other models of the Hall system as well. Finally, in Fig. 
4 we show the S 2 (7V) obtained from the last, band edge 
energy interval. Here, one has to modify the N such 
as to account for the finite size of the sample, that is, 
N —> Ni = N(l — N I K). Assuming as before a power 
law dependence for S 2 (#i), the best fit is obtained for 
a power of 7 = 0.97 ± 0.005 and a coefficient whose nat- 
ural logarithm is C = 0.007 ± 0.02. This is extremely 
close to the prediction of the Poisson ensemble but not 
identical, indicating that even in the tails of the Landau 
band there are a few states with relatively large £. 

The finite size correction used for the data of Fig. 4 
assumes that the levels are uncorrelated. In fact, the 
effect of this correction is not too large. Specifically, 
the best power law fit to Y, 2 (N) is 7 = 0.953 ± 0.005 
and C = 0.03 ± 0.02. This however, is because in Fig. 
4, N only goes up to about 60. On the other hand, 
the analogous finite size correction might significantly 
modify Fig. 3 where N reaches almost 350. In particular, 
it could modify the value of 71 . Since it is known that 
levels in this range are strongly correlated, it is not clear 
what the proper finite size correction should be. It would 
be interesting though to study the dependence of 71 on 
the system size. 
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3. Conclusions 

In summary, we have shown that the predictions of Refs. 
1-2 hold for the spectral statistics at the center of the 
lowest Landau band to the accuracy of our numerics. 
Moreover, it was pointed out that aside from the small 
and large s behavior of P(s), there is no real understand- 
ing of the spacings distribution interpolating between the 
metallic and the insulating regimes. In fact, the three 
most natural theoretical distributions, Pv(s), Ppm(s) 
and Pr(s), were shown to disagree with the numerical 
P(s). It should be instructive to further compare with 
the distribution derived by Haake and Lenz 21 from the 
way in which the Dyson gas relaxes to equilibrium from 
the configuration corresponding to a Poisson ensemble. 
However, in all these distributions one makes no explicit 
use of the fact that the eigenvectors of the system are 
localized. The simplest way to incorporate the presence 
of localization is in the spirit of renormalization theory. 
Namely, one assumes that the system is composed of in- 
dependent blocks of size ^ for which the corresponding 
spectrum has a Wigner P(s). The spectrum of the entire 
system is then a random superposition of the block spec- 
tra and the P(s) of such superpositions can be derived. 22 
Unfortunately, it turns out that the resulting P(s) has 
a nonvanishing value at s = in contradiction with the 
observed form. Recently, a new ensemble which has a 
preferred basis and thus, is reminiscent of the situation 
occurring in the case of localization, has been shown to 
be analytically tractable. 23 However, the P(s) for this 
ensemble is not yet known. 

Finally, we have shown that for very large N, £ 2 (A) 
grows with an unexpected power that is significantly 
larger than unity. 
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Fig. 1. The spacings distribution for the Nb = 119 levels 
closest to the center of the Landau band (•). The line is the 
corresponding best fit in the entire < s < 4 range of a) 
Pv(s), b) the Pandey-Mehta distribution, and c) the Robnik 
distribution. 
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Fig. 2. The spacings distribution for the Nb = 119 levels 
closest to E = 2 which, in turn, is close to the edge of the 
band (•). The best fitting Robnik distribution (solid) and 
the Poisson distribution (dashed) are also shown. 
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Fig. 3. The number variance, E 2 (7V), for the central en- 
ergy interval, \E\ < 0.5 (+). The solid lines are the corre- 
sponding power law fits in the three different N ranges (see 
text). 
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Fig. 4. Same as in Fig. 3 for the energy interval, 
1.5 < E < 2.5. Here however, the finite size corrected av- 
erage number of levels, Ni, is used as variable, instead of 
N. 



Table 1. The values of a and the corresponding \ 2 f° r 
gradually decreasing energy intervals around the center of 
the Landau band in which the P(s) was computed. 



N B 


a 


x 2 


119 


0.30 ± 0.02 


13.9 


59 


0.40 ± 0.03 


19.5 


29 


0.43 ± 0.05 


21.2 


19 


0.37 ± 0.06 


17.3 


9 


0.30 ± 0.08 


13.3 
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